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Verfahren zur Erfassung zeitabhangiger Moden 
dynamischer Systeme 



Zusammenf assung 



Bei einem Verfahren zur Erfassung der Moden eines dynamischen 
Systems mit einer Vielzahl von Moden, die jeweils einen Satz 
a(t) charakteristischer Systemparameter besitzen, wird eine 
Zeitreihe mindestens einer Systemvariablen x(t) einer 
Modellier.ung wie z. B. einer Schaltsegmentierung unterzogen, 
die dazu eingerichtet ist, in jedem Zeitabschnitt einer 
^^:gpyorbestiinmten Mindestlange fur jede Systemvariable x(t) ein 
vorbestimmtes Vorhersagemodell wie z. B. ein neuronales 
Netzwerk fUr eine entsprechende Systemmode zu erfassen, wobei 
nach der Modellierung der Zeitreihe eine Drif tsegmentierung 
erfolgt, bei der in jedem Zeitabschnitt, in dem das System von 
einer ersten Systemmode zu einer zweiten Systemmode Ubergeht, 
^ine Folge von gemischten Vorhersagemodeilen erfafit wird, die 
durchr eine lineare, paarweise Oberlagerung der 
Vorhersagemodelle der zwei Systemmoden gegeben ist. 
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Verfahren zur Erfassung zeitabhangiger Moden 
dynamischer Systeme 

Beschreibung 

Die Erfindung betrifft ein Verfahren zur Erfassung dynamischer 
Systeme, die durch im Zeitverlauf nichtstationare 
Systemparameter charakterisierbar sind, insbesondere ein 
Verfahren zur Segmentierung von Zeitreihen von Mefigrolien 
(Variablen) dynamischer Systeme und zur Identif izierung der 
die Segmente charakterisierenden Systemparameter (Moden) . 

Unter einem dynamischen System wird hier insbesondere jede 
Erscheinung verstanden, deren Zeitverlauf durch eine diskrete 
Abb il dung vom Typ 

x(t+l) = f^^,, (a(t) ) (0.1) 

darstellbar ist. Es werden jedoch auch Systeme mit mehreren 
(z. B. zwei) simultan auf genommenen Zeitreihen x, y gemaJi 

y(t+x) = f^^,, (x(t) ) (0.2) 

betrachtet. Dabei bezeichnen a(t) einen Satz 

charakteristischer Systemparameter, x einen Zustand, der im 
allgemeinen einen Vektor in einem mehrdimensionalen 
Zustandsraum bildet, und y einen zeitlich verschobenen 
Zustand. Der Zustandsraum wird von Variablen aufgespannt, die 
z^. B. physikalische, chemische, biologische, medizinische, 
geologische, geometrische, numerische und/oder 
prozeiitechnische Groiien sein konnen. 

Die Anzahl der Systemvariablen, die zusammen mit der Dynamik f 
das System beschreiben, entspricht der Dimension des 



Zustandsraumes, Hier werden Systeme betrachtet, deren 
Parameter a ebenfalls zeitlich veranderlich sein konnen. Ein 
gegebenes System mit zeitlich unveranderlichen Parametern a 
wird im folgenden auch als Mode bezeichnet. 

Beobachtbare oder meBbare Systemvariablen (MefigroBen) bilden 
erfaBbare Zeitreihen oder Datenstrome, die fur die jeweilige 
Abfolge von Systemmoden charakteristisch sind. Sind die 
Systemparameter innerhalb der 'Zeitreihen iiber bestimmte 
Zeitabschnitte (Segmente) unveranderlich, so kann die 
Zeitreihe entsprechend den jeweils bestehenden Systemmoden 
unterteilt werden (Segmentierung) und jedes Segment einer 
Systemmode zugeordnet werden ( Identif izierung) . 

Viele Vorgange in der Natur wie auch bei technischen 
Anwendungen wtirden sich dann vorhersagen und/oder 
kontrollieren lassen, wenn die ihnen zugrundeliegenden 
dynamischen Prozesse mathematisch modelliert werden konnen. 
Die Analyse und Charakterisierung praktisch gegebener 
dynamischer Systeme wird oft dadurch erschwert, dafi sich die 
Systemmoden wahrend der Beobachtung andern, Beispiele hier fur 
sind allmahliche Veranderungen, die sich in Drifts oder Trends 
der Systemparameter widerspiegeln, oder auch spontane oder 
abrupte Veranderungen in der Dynamik komplexer Systeme, . etwa 
.^■enn sich Konf igurationen spontan oder von auBen getrieben 
plotzlich andern, 

Ein Beispiel fur ein betrachtetes System ist die Erzeugung von 

Sprachsignalen im Mund-Rachen-Raum, bei der das System standig 

seine Konf iguration und somit seine Mode andert. Es besteht 

ein starkes Interesse daran, die Moden zu erfassen und zu 
•It 

identif izieren, die einem beobachteten Zeitverlauf einer 
Variablen (im Beispiel: Luf tdruckschwankungen) zugrundeliegen, 
auch xm verbesserte Vorhersagen liber das betrachtete System 
treffen oder dieses besser kontrollieren zu konnen. 



Prinzipiell konnen dynamische Systeme anhand voia gemessenen 
Signalen analysiert werden und es sind Reihe von Methoden 
bekannt, aus Zeitreihen Modelle zu gewinnen, die sich fur die 
Vorhersage und eine Kontrolle des Systemverhaltens eignen. So 
ist bekannt, daU sich der Zustand eines dynamischen Systems 
durch Erfassung der Zeitabhangigkeit beobachtbarer Mefigroiien 
modellieren lafit. Diese Modellierung erfolgt gema5 einem 
ersten Ansatz durch eine Rekonstruktion des Zustandsraumes 
mittels sogenannter Zeit-Verzo^gerungs-Koordinaten, wie es z. 
B. von N. H. Packard et al. in "Physical Review Letters" (Bd. 
45, 1980, S. 712 ff.) beschrieben ist. Auf der Basis einer 
solchen Rekonstruktion laflt sich dann nur ein einziges 
(globales) Modell f fur die Dynamik finden. Die globale 
Rekonstruktion des Systems ist ferner nachteilig, da bei 
Anwendungen fur vieldimensionale Systeme eine Vielzahl von 
Eingangsgrolien als Randbedingungen vorab bekannt sein mussen 
und/oder aufgrund der hohen Dimensionalitat das System sich 
praktisch nicht mehr schatzen (erfassen, abbilden) laBt 
und/oder ein ubermafiig hoher, nicht praktikabler Rechenaufwand 
entsteht • 

AuBerdem ist dieses Verfahren fur den Fall zeitlich 
veranderlicher Parameter allgemein nicht anwendbar. Die 
Analyse und Modellierung von dynamischen Signalen wird jedoch 
haufig dadurch erschwert, daiJ die zugrundeliegenden Systeme 
Ivlj^^^Sich in der Zeit in wesentlichen Parametern verandern. 

Beispiele sind z. B. Signale aus der Medizin, bei- denen ein 
Organ wie das Herz oder das Gehirn viele dynamische Moden hat, 
die sich abwechseln, oder etwa den Sprachsignalen, bei denen 
das erzeugende System, .der Mund-Rachen-Raum, in der Zeit 
of f ensichtlich verschiedene Konf igurationen durchlauft. 

Gemafi einem weiteren Ansatz ist daher aus der Publikation von 
K, Pawelzik, J. Kohlmorgen und K.-R. Muller in "Neural 
Computation" (Bd. 8, 1996, S. 340 ff.) bekannt, Datenstrome 
entsprechend zunachst unbekannten, zeitlich wechselnden 



Systemmoden durch Simulierung mit mehreren, untereinander in 
Konkurrenz stehenden Modellen zu segmentieren. Die Modelle 
werden bevorzugt durch neuronale Netzwerke gebildet, die 
jeweils fUr eine Dynamik charakteristisch sind und nach 
vorbestimmten Trainingsregeln urn die Beschreibung der 
einzelnen Punkte des Datenstroms konkurrieren . 

Mit dieser Methode gelingt es, eine Zeitreihe in Abschnitte 
quasistationarer Dynamik zu zerlegen und simultan Modelle fur 
diese System-Moden aus den Zeitreihen zu identif izieren. 

Die Segmentierung gemaij K. Pawelzik et al., zu der unten 
Einzelheiten angegeben werden, erlaubt die Zuordnung von 
Segmenten zu bestimmten Systemdynamiken oder -moden und fuhrt 
zu einer Erfassung des Datenstroms als Vorgang mit einem 
diskreten Schalten (sog. "switching") zwischen den Moden. 
Diese Beschreibung der Parameterdynamik komplexer Systeme 
stellt zwar gegenuber der o. a. globalen Modellierung 
hinsichtlich der Genauigkeit und der Segmentierung 
verschiedener Systemzustande einen Fortschritt dar. Allerdings 
kann der Obergang zwischen verschiedenen Systemzustanden nicht 
ausreichend beschrieben werden. Es hat sich insbesondere bei 
der Analyse realer Systeme z, B. bei medizinischen Anwendungen 
gezeigt, daB die Segmentierung auf bestimmte Falle mit 
;noglichst klaren Modenunterschieden und geringem Rauschen 
^l^^eschrankt, allgemein jedoch bei zeitlichen Veranderungen der 
erzeugenden Systeme unzuverlassig ist. 

Solche zeitlichen Veranderungen der erzeugenden Systeme machen 
die beobachtbaren Signale nichtstationar und fUhren dazu, dafl 
sich die Systeme im allgemeinen nicht mehr durch einheitliche 
Medelle beschreiben lassen. Erfolgen diese Veranderungen der 
Systeme plotzlich, so spricht man von Sprungprozessen. 

Die Aufgabe der Erfindung ist es, verbesserte Verfahren zur 
Erfassung der Moden dynamischer Systeme mit nichtstationaren 



Systemparametern anzugeben, mit denen die Beschrankungen 
herkommlicher Verfahren uberwunden werden konnen und die es 
insbesondere ermoglichen, automatisch mit praktikablem 
Bearbeitungsaufwand und hoher Zuverlassigkeit Segmentierungen 
und Identifizierungen von Zeitreihen mit einer erhohten Zahl 
von Einzelheiten vorzunehmen. 

Diese Aufgabe wird durch das Verfahren mit den Merkmalen gemali 
den Patentanspruch 1 .gelost. Vlarteilhaf te Ausfuhrungsformen 
der Erfindung ergeben sich aus den abhangigen Anspruchen. 

Die Erfindung basiert auf der Idee, Obergange zwischen 
verschiedenen Moden eines dynamischen Systems als 
Zwischenmoden des Systems aufzufassen, die paarweise lineare 
Interpolationen der Ausgangs- bzw. Endmoden des Ubergangs 
darstellen. Die betrachteten dynamischen Systeme gehen eher 
allmahlich von einer Mode in die andere uber, statt abrupt 
zwischen Moden zu schalten. Die Erfindung zielt darauf ab, 
solche Ubergange zwischen dynamischen Moden in Signalen und 
die Moden zu identif izieren . 

Es wird daher bei einem Verfahren zur Erfassung der Moden 
dynamischer Systeme z. B. nach einer Schaltsegmentierung einer 
Zeitreihe mindestens einer der Systemvariablen x(t) des 
_ . eine Drif tsegmentierung vorgenommen, bei der in jedem 

'^%|i|^eitabschnitt/ in dem das System von einer ersten Systemmode s^ 
zu einer zweiten Systemmode s^ ubergeht, eine Folge von 
gemischten Vorhersagemodellen g^ erfaiit wird, die durch eine 
lineare, paarweise Uberlagerung der Vorhersagemodelle f^^^ der 
zwei Systemmoden s^^ ^ gegeben ist . 

G*egenstand der Erfindung ist auch eine Vorrichtung zur 
Erfassung eines dynamischen Systems mit einer Vielzahl von 
Moden Si mit jeweils charakteristischen Systemparametern a(t). 
Die Vorrichtung enthalt eine Einrichtung zur Aufnahme einer 
Zeitreihe mindestens einer der Systemvariablen x(t) des 



Systems, eine Schaltsegmentierungseinrichtung, die dazu 
eingerichtet ist, in jedem Zeitabschnitt einer vorbestiinititen 
Mindestlange fUr die Systemvariablen x{t) ein vorbestimmtes 
Vorhersagemodell fur eine entsprechende Systemmode s^ zu 
erfassen, und eine Drif tsegmentieriingseinrichtung, mit der in 
jedem Zeitabschnitt, in dem das System von einer ersten 
Systemmode s^ zu einer zweiten Systemmode s^ Obergeht, eine 
Folge von gemischten Vorhersagemodellen g^ erfafit wird. Die 
erf indungsgemafle Vorrichtung lcann_ f erner Einrichtungen zur 
Einstellung von Interpolations- und Segmentierungsparametern, 
Vergleicherschaltungen zur Verarbeitung von Vorhersagef ehlern 
von Vorhersagemodellen, Anzeige- und Signalisierungsein- 
richtungen und Speichereinrichtungen enthalten. Die 
erf indungsgemafle Vorrichtung kann ein Monitor fUr 
physiologische Daten oder physikalische oder chemische 
Prozeiiparameter sein, 

Mit der Erfindung steht ein Instrument zur Verfiigung, welches 
ein grofies Anwendungspotential in vielen medizinischen, 
wissenschaf tlichen und technischen Bereichen hat, Mit der 
Segmentierung von Signalen und einer entsprechenden 
Identif izierung der zugrundeliegenden Dynamik eroffnen sich 
neue Moglichkeiten der Vorhersage und der Kontrolle auch 
wesentlich nichtstationarer Systeme. 

^2\nwendungen der Erfindung haben gezeigt, daB sich 
kontinuierliche Obergange zwischen Systemmoden- sicher 
identif izieren lieUen und daB die zugrundeliegenden Dynamiken 
durch die Modelle mit einer Prazision beschreibbar sind, die 
in vielen Fallen eine Vorhersage des Systemverhaltens 
ermoglicht. Mit der Erfindung lassen sich in vielen Fallen 
nichtstationarer Prozesse Modelle identif izieren, die sich 
auch fur eine Kontrolle der Prozesse eignen, die ohne eine 
Beriicksichtigung der Nichtstationaritat nicht moglich waren. 



Ausfuhrungsformen und weitere Vorteile der Erfindung warden im 
folgenden unter Bezug auf die beigefiigten Zeichnungen 
beschrieben. Es zeigen: 

Fig. 1: Kurvendarstellungen zur Illustration eines ersten 
Segmentierungsschrittes des erf indungsgemaiien Verfahrens; 

Fig. 2: Kurvendarstellungen zur Illustration eines weiteren 
Segmentierungsschrittes des er.f.indungsgemaiien Verfahrens 

Fig. 3: Kurvendarstellungen zur Segmentierung von 
Blutregulierungsdaten nach dein erf indungsgemaiien Verfahren; 
und 

Fig. 4: Kurvendarstellungen zur Segmentierung von EEG-Daten 
nach dem erf indungsgemaiien Verfahren. 

Im folgenden werden zunachst Einzelheiten der Erfindung unter 
Bezug auf die Figuren 1 und 2 und anschlieliend praktische 
Anwendungsbeispiele erlautert, Es ist dem Fachmann 
ersichtlich, daii die Erfindung nicht auf die 

Anwendungsbeispiele beschrankt ist, sondern entsprechend auch 
in anderen Gebieten angewendet werden kann, wie sie 
beispielhaft welter unten aufgefuhrt sind. 

^^^gp^'(l) Erfassung von Drif tubergangen in nichtstationaren 
Zeitreihen 

Erf indungsgemaii werden nichtstationare Zeitreihen mit einer 
zweischrittigen Prozedur erfaiit, bei der zunachst eine 
geeignete Modellierung und anschlieliend eine sogenannte 
Drif tsegmentierung erfolgt. Die Modellierung ist dazu 
eingerichtet , in jedem Zeitabschnitt einer vorbestimmten 
Mindestlange fur jeden Systemparameter ein vorbestimmtes 
Vorhersagemodell fur eine entsprechende Systemmode zu 
erfassen. Hierzu erfolgt vorzugsweise eine herkommliche 
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Schaltsegmentierung, wie sie beispielsweise aus der 
Publikation von K. Pawelzik et al . in "Neural Computation" 
(Bd. 8, 1996, S. 340 ff.) bekannt ist. Die Modellierung kann 
aber auch durch eine andere, in Bezug auf die gewonnene 
Systeminf ormation zur Schaltsegmentierung aquivalente Prozedur 
erfolgen, die an eine konkrete Anwendung z. B. bei bekannten 
reinen Moden oder Randbedingungen angepalit ist. 

Die Schritte der Schalt- und Dx.if tsegmentierung werden im 
folgenden naher erlautert. Hinsichtlich der 

Schaltsegmentierung wird hiermit der Inhalt der Publikation 
von K. Pawelzik et al. ausdrUcklich vollstandig in die 
vorliegende Beschreibung einbezogen . 

(i) Schritt 1 (Schaltsegmentierung) 

Die Schaltsegmentierung dient der Ermittlung 
charakteristischer Prediktoren, die zur Beschreibung der 
Systemmoden geeignet sind. Die Schaltsegmentierung kann 
entweder an einer Trainings-Zeitreihe oder an der zu 
untersuchenden Zeitreihe durchgefuhrt werden. In beiden Fallen 
konnen die ermittelten Vorhersagemodelle oder Prediktoren fur 
weitere, unbekannte Zeitreihen verwendet werden. 

J.s wird ein dynamisches System mit einer endlichen Zahl N 
"^^i^>??^nterschiedlicher Moden betrachtet. Fur die j-te Mode ist ein 
Wert (allg.: Vektor oder Satz) CL^{t) eines beobachtbaren 
Systemparameters charakteristisch, der jeweils mit einer 
Funktion f^^^, (i = 1,...,N) aus einem Satz von N Funktionen f 
modelliert werden soil'. Es wird die Zeitreihe {x^} = (t) der 
Systemvariablen betrachtet und zeitabhangig jeweils nach der 
Funktion f^^^, gesucht, fur die {yj = y^(t) = f,^^^(x^(t)) eine 
neue Zeitreihe vorherzusagender Punkte y^it) darstellt, die in 
Bezug auf die Systemmoden qualitativ die selben Eigenschaf ten 
hat wie {x^} . Durch den Wechsel der Modellfunktion f im 
Zeitverlauf wird die Schaltsegmentierung gefunden, die die 




Zeitreihe {x^} entsprechend der wechselnden Systemmoden 
unterteilt . 

Die Funktionen f werden als Prediktoren (oder: 
Vorhersagemodelle, Expertenfunktionen) aus einem Satz von 
Netzwerken mit veranderlichen Parametern durch ein geeignetes 
Trainingsprogramm abgeleitet, bei dem simultan sowohl die 
Parameter der Netzwerke als auch die Segmentierung ermittelt 
werden. Der Begriff "Netzwerk". wird hier fur alle moglichen 
geeigneten Modellfunktionen verwendet, also vorzugsweise fiir 
neuronale Netzwerke, aber auch z. B. fur Polynome oder lineare 
Funktionsapproximationen, Die optimale Wahl eines neuronalen 
Netzwerkes erfolgt in Abhangigkeit von der spezif ischen 
Anwendung. Es werden vorzugsweise Netzwerke mit einer 
schnellen Lernf ahigkeit , wie z. B, sogenannte RBF-Netzwerke 
(Radial Basis Function Network) vom Moody-Darken-Typ 
verwendet . 

Das Training erfolgt unter der Voraussetzung, daB die 
Systemmoden nicht mit jedem Zeitschritt wechseln, sondern eine 
geringere Schaltrate aufweisen, so daii eine Systemmode tiber 
mehrere Zeitschritte erhalten bleibt. Die angenommene Grenze 
der Schaltrate bzw. Anzahl der Zeitschritte, iiber die eine 
Systemmode erhalten bleibt, ist zunachst ein freier 
. "ingangsparameter und kann anwendungsabhangig in geeigneter 
^i^SS^eise beispielsweise in Abhangigkeit vorgegebener 

Erf ahrungswerte oder nach einer Parameteranpassungsstrategie 
ausgewahlt werden. Bei der Parameteranpassungsstrategie kann 
vorgesehen sein, fiir die Schaltrate einen Ausgangswert 
vorzugeben und mit diesem einen Vorhersagef ehler (siehe unten) 
zu ermitteln. Falls die Schaltrate zu hoch oder zu niedrig 
gewahlt ist, fiihrt eine entsprechende Ober- oder 
Unterspezialisierung zu einem zu hohen Vorhersagef ehler , Im 
weiteren Verlauf der Anpassung kann dann die Schaltrate 
optimiert werden, bis der mittlere Vorhersagef ehler unterhalb 
vorbestimmter Grenzen liegt. 
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Das Training erfolgt durch Maximierung der Wahrscheinlichkeit 
W, dafl der Satz von Netzwerken die Zeitreihe {x^} erzeugen 
wiirde. Es handelt sich um ein Training mit konkurrierendem 
Lernen (engl.: "competitive learning"), wie es im einzelnen in 
der Abhandlung "Introduction to the theory of neural 
computation" von J. Hertz et al* (Addison-Wesley Publishing " 
Company 1991, insb. Kap. 9: "Unsupervised competitive 
learning") beschrieben wird. Aus dieser Abhandlung ist die 
anwendungsabhangige Implementierung eines solchen Trainings 
ableitbar. Die Trainingsregel des konkurrierenden Lernens auf 
der Grundlage des beim Lernen auftretenden Fehlers ist gemaB 
(1) darstellbar. 



Diese Trainingsregel stellt sicher, dafl die 
Lerngeschwindigkeit (Verbesserung der Parameter) fur die 
Funktionen f mit geringstem Abstand zum Zielwert y am hochsten 
ist . 

Figur 1 zeigt das Ergebnis der Schaltsegmentierung am Beispiel 
der Analyse einer chaotischen Zeitreihe {x^} mit x^.^^ = f (x^) 
zwischen den vier Moden: 



fi{x) = 4x(l-x) fiir X G [0, 1] - 
f2(x) = f,(f,(x) ) 

f3(x) = 2x ftir X e [0, 0.5) bzw, 

f3{x) = 2(l-x) fur X e* [0.5, 1] 
f,(x) = f3(f3(x) ) 

Fur die ersten 50 Zeitschritte wird mit einem Startwert 
Xq = 0.5289 zunachst f^ angewendet. AnschlieBend erfolgt ein 
Obergang (Einzelheiten siehe (ii) ) zur Mode f2, die nach 
Schritt 100 bis Schritt 150 stationar wird. Entsprechend wird 



dlog W 



OC 




(y-fi) 



(1) 



ab Schritt 200 bzw. Schritt 300 jeweils fur 50 Schritte die 
Mode f3 bzw. eingenommen. Danach erfolgt der Riickubergang 
auf f^. Figur la zeigt einen Ausschnitt (Schritte 300 bis 450) 
des Zeitverlaufs der Zeitreihe {x^} mit x^^^ = f(x^.). 

Die Segiuentierung der ersten 450 Zeitschritte mit 6 
Prediktoren t^, i = 1,...,6 (RBF-Netzwerke vom 
Moody-Darken-Typ) ist in Figur lb dargestellt. Das Training 
ergibt eine Spezialisierung vojx. vier der Prediktoren (6, 2, A, 
3) jeweils auf eine der o. a. vier Moden. Die stationaren 
Bereiche liegen bei den Intervallen [0, 50] und [400, 450] 
(fj, [100, 150] (f^), [200, 250] (f3) und [300, 350] (fj. Die 
tibrigen zwei Prediktoren (3, 5) tiaben sich auf die 
Obergangsbereiche zwischen den Moden spezialisiert . Dies zeigt 
den Naciiteil der herkommliclien Sclialtsegmentierung, bei der im 
Fall von Obergangen der eritsprechende Zeitbereicli ohne 
adaquate Besclireibung mehrfacli unterteilt wird. 

Anstelle des hier beschriebenen sog. "Hard competition"- 
Trainings, bei dem bei einem Trainingssciriritt jeweils nur ein 
Vortiersagemodell optimiert wird ("winner takes all") kann auch 
vorgeselien sein, den Grad der Konkurrenz im Ralimen eines sog. 
"Soft competition"-Trainings zu verandern, wie es im einzelnen 
in der Publikation von K. Pawelzik et al. besclirieben ist. 

(ii) Scliritt 2 (Drif tsegmentierung) 

Beim zweiten Schritt werden die Obergange (sog. Drif ten, 
nicht-abruptes Obergehen, gleitendes Wechseln) zwischen den 
Systemmoden berlicksichtigt . Mit der Erfindung wurde als 
wichtige Voraussetzung fur die Drif tsegmentierung 
herausgefunden, daiJ der Obergang von einer ersten Systemmode 
direkt in eine zweite Systemmode und nicht uber eine dritte 
Systemmode erfolgt. Das Drif ten zwischen Systemmoden wird 
somit wie folgt durch eine Oberlagerung von (oder paarweise 



lineare Interpolation zwischen) genau zwei Moden modelliert . 
Dabei treten gemischte, ggf, abgestufte Zwischenmoden auf, die 
jedoch nicht eigene (reine) Systemmoden sind. 

Es wird ein Satz von P reinen Systemmoden, die jeweils durch 
ein Netzwerk k(s), s e reprasentiert werden, und ein Satz 
von M gemischten Systemmoden betrachtet, die jeweils durch 
eine lineare Oberlagerung von zwei Netzwerken i(s) und j (s) , 
s € M, reprasentiert werden. Das Modellnetzwerk g^ fur eine 

gegebene Mode seS, S = PoM wird gemafi Gleichung (2) 
dargestellt . 




In (2) ist X der Vektor (x^, x^^, . . , , x^.^^.^,^) der 
Zeitverzogerungskoordinaten der Zeitreihe {x^} und f^^^ sind 
Prediktoren, die gemali der o. a. Schaltsegmentierung ermittelt 
wurden. Dabei ist m eine Einbettungsdimension und x der 
Verzogerungsparameter der Einbettung. Die Einbettungsdimension 
ist die Dimension des Phasenraumes, in dem das System 
betrachtet wird und in dem die Modelle operieren. 



V ^ur jede gemischte Systemmode sind zwei Parameter a, b 
zusammen mit zwei Netzwerkindizes i, j charakteristisch . Zur 
Vereinf achung des Rechenaufwandes wird die Zahl'der gemischten 
Moden begrenzt. Mit 0 < a(s) < 1 und b(s) = 1 - a(s) wird 
eine endliche Zahl von Werten a{s) definiert. Zur weiteren 
Vereinf achung werden gleiche Abstande zwischen den Werten a(s) 
gemafl Gleichung (3) gewahlt. 

ar = R I ^ niit r = 1,..., R (3) 



R entspricht der Anzahl der zugelassenen Zwischenmoden und 
wird auch als Auflosung oder Abstufung der Interpolation 
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zwischen den reinen Moden bezeichnet. Die Auflosung R kann 
einen beliebigen Wert annehmen, wird jedoch anwendungsabhangig 
zur Erzielung einer optimalen Systembeschreibung (insbesondere 
bei stark verrauschten Vorgangen) und praktikabler 
Rechenzeiten insbesondere unter Berucksichtigung der oben 
genannten Schaltrate ausreichend niedrig gewahlt. Bei 
praktischen Anwendungen (siehe unten) kann vorgesehen sein, 
daB die Auflosung R manuell von einem Bediener oder 
automatisch von einem Stellkrei-S in Abhangigkeit von einem 
vorliegenden Analyseergebnis und einem Vergleich mit einem 
vorgegebenen Schwellwert gewahlt werden. 

Die Gesamtzahl der gemischten Moden betragt bei einer 
gegebenen Auflosung R zwischen jeweils zwei Netzen 
|M| = R • N • (N-l)/2. Bei dem o, a. Beispiel betragt somit bei 
N = 8 reinen Moden und einer Auflosung R = 32 und Gesamtzahl 
der gemischten Moden |M| = 896. Zur Bestiramung der Gesamtzahl 
der Systemmoden kommen noch die 8 reinen Moden hinzu. 

Die Drif tsegmentierung umfaiit nun die Suche nach einer 
Segmentierung mit den reinen und gemischten Systemmoden (a, b, 
R) , die in Bezug auf den Vorhersagef ehler der Moden der 
gesamten Zeitreihe optimiert ist. Die Prediktoren werden so 
ausgewahlt, daB jedem Element der Zeitreihe einer der Moden 
';.us der Gesamtzahl der Systemmoden zugeordnet werden kann. Der 
Sv^'orhersagef ehler ist die Abweichung einer 
Prediktorenvorhersage vom tatsachlichen Element* der zu 
untersuchenden Zeitreihe. Fiir die zu untersuchende Zeitreihe, 
die nicht mehr unbedingt die Trainings-Zeitreihe ist, mit der 
bei der Schaltsegmentierung die angepaliten Netzwerke oder 
Prediktoren ermittelt wurden, wird fiir jeden Zeitschritt mit 
jisdem der Prediktoren eine Vorhersage ermittelt, woraus sich 
eine zeitabhangige Matrix der Prediktorenvorhersagen ergibt, 
aus der ein mittlerer Vorhersagef ehler fiir willkiirlich 
gewahlte Segmentierungen ableitbar ist. Die Segmentierungen 
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mit deiti geringsten Vorhersagef ehler ist die gesuchte 
Drif tsegmentierung. 

Die Suche nach der Segmentierung mit dem geringsten 
Vorhersagef ehler kann mit jeder geeigneten Such- oder 
Iterationstechnik erfolgen, Vorzugsweise wird eine dynamische 
Programmiertechnik gewahlt, die aquivalent zum 

Viterbi-Algorithmus fUr HM-Modelle (sog. Hidden Markov Models) 
ist. Einzelheiten hierzu sind jD-eispielsweise in der 
Publikation "A Tutorial on Hidden Markov Models and Selected 
Applications in Speech Recognition" von L. R, Rabiner in 
"Readings in Speech Recognition" (Hrsg. A. Waibel et al., San 
Mateo^ Morgan Kaufmann, 19*90, S. 267-296) beschrieben. Die 
Drif tsegmentierung ist im Rahmen der HM-Modelle die 
wahrscheinlichste Modensequenz, die die zu untersuchende 
Zeitreihe erzeugt haben kdnnte . Als Nebenbedingung wird dabei 
die Moglichkeit der Moden-Anderungen Ober die T-Funktion 
(siehe unten) eingeschrankt . 

Das Ziel der Anpassung ist die Angabe einer optimalen Sequenz 
von Netzwerken bzw, linearen Mischungen aus diesen. Eine 
Sequenz ist optimal, wenn die sogenannte Energie- oder 
Cost-Funktion C* der Vorhersage minimiert ist. Die 
Cost-Funktion C* setzt sich aus der Summe der quadratischen 
"^ehler der Vorhersage und den Cost-Funktionen der 
fe^;^Siodenlibergange der Sequenz zusammen. Die Ableitung der 

Cost-Funktion C* jeweils zwischen zwei Zeitpunkten tg und t^^ 
erfolgt induktiv, indem zunachst von einer Start-Cost-Funktion 
gemaJi Gleichung (4) ausgegangen wird. 



CJto) = 8, (to) 



(4) 



wcbei 



e,(t) = (X, - g,(x,.J)^ 



(5) 



der quadratischen Fehler der Vorhersage der reinen oder 
gemischten Moden g ist. 



Fiir den Induktionsschritt von t - 1 auf t wird die 
Cost-Funktion gemaii Gleichung (6) fur alle s e S berechnet: 

CJt) = ejt) + xnin {Cs{t-1) + T(s, s)}, t=to+l , . . . tmax (10) 

Dabei ist T(s, s) die Cost-Funktion des tJbergangs von einer 
Mode S zu einer Mode s. 

Die optimale (minimale) Cost-Funktion C* ist dann: 

C* = min {Cs (t^^J } (11) 

Die Funktion T entspricht bei den HM-Modellen den 
Ubergangswahrscheinlichkeiten und kann anwendungsabhangig 
geeignet gewahlt werden. Es kann beispielsweise vorgesehen 
sein, abrupte Schaltiibergange und gleitendes Driften zwischen 

zwei Netzen zuzulassen und alle anderen Obergange durch T = oo 
auszuschlielien. 

Die Drif tsegmentierung ergibt sich aus der ermittelten 
optimalen Sequenz von Netzwerken bzw. linearen Mischungen aus 
diesen, indem die Moden, die C* ergeben, zuruckverf olgt und in 
Abhangigkeit von der Zeit erfalit werden. 

'^§0±m AnschluB an die Drif tsegmentierung kann sich ein 
zusatzlicher Schritt der Reduzierung der Zahl der ■ zur 
Modellierung verwendeten Netzwerke anschlieiien, der im 
einzelnen unten beschrieben wird. 

Abschlieiiend findet eine Identif izierung der segmentierten 
Moden statt, indem jedem Prediktor oder Vorhersagemodell die 
zugehorige Systemmode zugeordnet wird. Diese Identif izierung 
erf olgt anwendungsabhangig . 
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Das Ergebnis der Drif tsegmentierung im Fall der oben unter 
Bezug auf Figur 1 erlauterten chaotischen Zeitreihe {xj mit 
vier Moden wird im folgenden unter Bezug auf Figur 2 
beschrieben. Die Drif tsegmentierung umfaiJt die Suche nach 
einem Verlauf a(t), der einen speziellen Pfad zwischen den 
reinen Moden liefert, fUr den der Vorhersagef ehler der 
gesamten Zeitreihe optimiert ist. 

Jewells nach den ersten 50 ZeLtschritte mit der Mode gemaB f^ 
erfolgt fUr 50 Schritte ein zeitlinearer Dbergang in die Mode 
gemaB t^. Der Obergang ist ein zeitabhangiges Drif ten gemaB den 
Gleichungen (12) : 



Entsprechende Obergange finden nach dem 150., 250. und 350. 
Schritt fur jeweils 50 Schritte statt. 

In Figur 2 ist die Besetzung der jeweiligen Moden entsprechend 
den ermittelten Netzwerken in Abhangigkeit von der Zeit 
(Zeitschritte [1200, 2400]) dargestellt. Aus 
Ubersichtlichkeitsgrunden sind die Obergangs- oder 
Orif tbereiche entsprechend ihren Zeitgrenzen und Ausgangs- 
^■0^^2M. Endmoden in Rahmen dargestellt, in denen der jeweilige 
Driftverlauf zwischen den Moden gepunktet ist. ._ Figur 2a zeigt 
bei einer Auflosung R = 32 (siehe Gleichung (3)) Obergange wie 
z. B. bei den Zeitschritten 1350 bis 1400 zwischen den Netzen 
2 und 4. Die Obergange sind wie gemaB Gleichung (12) erwartet 
linear. Bei einer geringeren Auflosung von R = 3 ergibt sich 
d4e in Figur 2b dargestellte Segmentierung. Abweichend vom 
linearen Driftverlauf sind die gepunkteten Obergange 
stuf enformig. Dennoch ist die Darstellung auch bei der 
geringeren Auflosung eine adequate Beschreibung des 



f(x,) = (1 - a(t))f,(x,) + a(t)f2(x,) 



mit 



(12) 



^(t) = t . - t\ = 50, t^ = 100 



dynamischen Verhaltens des Systems, wie der Vergleich der 
zeitlichen Lage der Moden und des Driftens zeigt. 

(2) Anwendungsbeispiele fiir die Erfassung von Drif tUbergangen 

(i) Blutzellenregulation im menschlichen korper 

Die Blutzellenregulation im menschlichen Korper stellt ein 
hochdimensionales chaotisches System dar, das durch die 
Mackey-Glass-Verzogerungs-Dif f erentialgleichung (13) (siehe 
auch in der o. a. Abhandlung von J. Hertz et al.) beschrieben 
werden kann. 



Erf indungsgemaii konnen Zeitreihen von physiologischen 
Parametern, die fur die Menge der roten Blutkorper 
charakteristisch sind, anwendungsabhangig segmentiert werden. 
Die Funktionsf ahigkeit der Segmentierung wird im folgenden 
beispielhaf t erlautert . 

Bei Vorgabe von zwei Moden A und B, die sich jeweils durch 
verschiedene Verzogerungsparameter t^ = 17 bzw. t^ = 23 
, j.nterscheiden, erfolgt mit einer Abtast-Zeitschrittgrolie 1=6 
'"'rw'nach 100 Schritten zunachst ein Obergang von A nach B. Der 
Obergang dauert 100 Schritte und ist eine Oberlagerung der 
Gleichung (13) mit den beiden Verzogerungsparameter t^ wahrend 
der Integration von Gleichung (13) . Die Oberlagerung wird mit 
einem exponentiellen Drif tparameter a (siehe Gleichung (2)) 
gemaii Gleichung (14) erzeugt. 



dx(t) 
dt 



= -O.lx(t) + 



0.2x(t-td) 
l+x{t-td} ^° 



(13) 




t = 



1 



100 



(14) 



In der Folge wiederholen sich alle 100 Schritte stationare 
Moden A oder B bzw. die jeweiligen Obergange. Fur jeden 
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RUck-Obergang nach einem Drif tObergang wird ein schaltartiges 
Umspringen angenoinmen • Figur 3a zeigt die entsprechende 
Zeitreihe iiber 300 Schritte. Die Drif tsegmentierung mit sechs 
Prediktoren auf der Grundlage von RBF-Netzwerken mit jeweils 
40 Basisfunktionen, einem Einbettungsparameter m = 6 und dem 
Verzogerungsparameter x = 1 (siehe Gleichung (2)) ergibt das 
in Figur 3b gezeigte Bild. Es zeigt sich die efwartete 
Segmentierung der Zeitreihe in stationare Moden und 
Drif tiibergange . - 



Allerdings zeigt sich, dafl sich jeweils 2 Netzwerke auf eine 
Mode spezialisiert haben (2, 3 => Mode A, 5, 6 => Mode B) , In 
siner solchen Situation kann es erf indungsgemaB vorgesehen 
sein, den zusatzlichen Schritt der Reduzierung der Zahl der 
zur Modellierung verwendeten Netzwerke vorzusehen. 

Der Reduzierungschritt umfaJit eine sequentielle Reduzierung 
der Netzwerkzahl, jeweils verbunden mit einer Bestimmung des 
mittleren Vorhersagef ehlers . Die Reduzierung (Abzug von 
Redundanznetzwerken) wird beendet, falls eine weitere 
Verringerung der Netzwerkzahl eine signifikante Erhohung des 
Vorhersagef ehlers bedeutet. Figur 3c zeigt das Ergebnis einer 
solchen Reduzierung. Der mittlere quadratische Fehler RMSE 
bleibt bei der Reduzierung urn ein, zwei, drei und vier 
■-..etzwerke jeweils konstant, wohingegen ein starker Anstieg fur 
den Fall auftritt, dafi nur noch mit einem Netzwerk modelliert 
wird. Dies bedeutet, dali das System optimal mit einer Zahl von 
Netzwerken modelliert wird, die gleich der Gesamtzahl von 
betrachteten Netzwerken, veirmindert um die Zahl von 
Redundanznetzwerken ist . 

Die adequaten Modellnetze werden dann erhalten, indem der 
RMSE-Wert ftir jede Netz-Kombination mit reduzierter Netzzahl 
berechnet wird. Die Netz-Kombination mit dem geringsten 
RMSE-Wert umfaflt die gesuchten Modellnet zwerke Oder 
Prediktoren- Figur 3d zeigt die Drif tsegmentierung nach dem 




Reduzierungsschritt . Demnach beschreiben die tibrig gebliebenen 
Prediktoren 2 und 5 das System vollstandig. 

(ii) Schlafdatenerf assung 

Eine weitere Anwendung der Erfindung liegt im Ber^ich der 
Analyse physiologischer Daten, die fiir den Ablauf von Schlaf- 
und Wachmoden von Lebewesen cbarakteristisch sind. Als 
Grundlage fiir sich anschlieiiende Prozeduren zur Erfassung von 
Schlaf storungen konnen Zeitreihen z. B. von EEG-Daten 
segment iert werden. 

Figur 4a zeigt im Vergleich die Ergebnisse einer herkommlichen 
Schaltsegmentierung (oben) , einer Drif tsegmentierung (Mitte) 
und einer "manuellen" Segmentierung (unten) eines 
medizinischen Fachmannes (Schlaf forscher) auf der Grundlage 
von Erf ahrungswerten am Beispiel eines Nachmittagsschlaf es 
eines gesunden Menschen. Die Schalt- und Drif tsegmentierungen 
erfolgen mit acht Netzwerken (netl . . . netS) an 
Einkanal-EEG-Daten x(t) (Figur 4b). In Figur 4a sind wie bei 
Figur 2 aus Obersichtlichkeitsgrunden Rahmen gezogen, die bei 
den Driftmoden verdeutlichen, zwischen welchen Netzen 
interpoliert wird. Die gepunktete Linie im Inneren der Rahmen 
. ;eigt den jeweils tatsachlichen Verlauf. Die manuelle 
;:i5egmentierung basiert auf der Beobachtung von physiologischen 
Signalen (z. B. EEG, EOG, EGG, Puis, Blutdruck; ' Atmung, 
Augenbewegungen) . Die Moden Wl, W2 bezeichnen zwei Wachmoden 
mit offenen bzw. geschlossenen Augen und die Moden SI, S2 
jeweils Schlaf zustande-. "n.a." und "art." beziehen sich auf 
nicht in Betracht gezogene Zustande bzw. Artefakte. 

Die Schaltsegmentierung zeigt ein vergleichsweise 
undif f erenziertes Bild, das nur grob mit den ubrigen 
Beobachtungen konsistent ist. So tritt beispielsweise ein 
Einschlafvorgang in alien drei Fallen bei t « 7000 auf. Die 



Driftsegmentierung ergibt jedoch mehrere Drif tUbergange, die 
zusatzliche Einzelheiten des Schlafverhaltens darstellen. Der 
"manuell" beobachtete Schlafbeginn zur Zeit t « 4000 wird 
durch einen exponentiellen Drif tubergang von Netz net7 
(Wachmoden-Prediktor) zum Netz net4 (Schlafmoden-Prediktor ) 
reprasentiert . Das Aufwachen setzt bei t « 9Q00 durch ein 
leichtes Rilckdriften ziim Netz net7 ein, das bis zum Erreichen 
des "manuell" ermittelten Aufwachpunktes t « 9500 gehalten 
wird. In dieser Situation erfolgt eine plotzliche Anderung des 
Wichtungsf aktors, so daii das Netz net7 eine groBere Wichtung 
erhalt. Nach t « 9800 (Augen of fen) besteht eine Mischung der 
beiden Wachmoden-Prediktoren net7 und net2. 

(iii) Weitere Anwendungen und Vorteile 

Figur 4a zeigt, daB mit dem erf indungsgemaiien Verfahren 
automatisch detaillierte Segmentierungen erzielt werden 
konnen, die bislang nur durch die Beobachtung komplexer 
Merkmalsbilder auf der Grundlage breiter Erfahrungen und 
Intuitionen zuganglich waren. Dieser Vorteil ist nicht nur in 
der Medizin, sondern auch in anderen Gebieten nutzbar, in 
denen groiie Datenmengen bei der Beschreibung komplexer 
iynamischer Systeme anf alien. Derartige Gebiete sind die 
■^l^l^^hysikalische, chemische und/oder biologische 

Verf ahrenstechnik, die Geologie, Meteorologie, - Klimatologie, 
die Spracherf assung u. dgl , . 

Erf indungsgemafie Verfahren haben die folgenden Vorteile. Das 
betrachtete System kann hochdimensional sein (10 oder mehr 
Dimensionen) . Die Erfindung erlaubt eine Reduzierung der 
Komplexitat eines solchen Systems durch Betrachtung von 
niedrigerdimensionalen Moden und von wechselnden Obergangen 
zwischen diesen, Der Einsatz von Vorhersagemodellen fur die 
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Segmentierung ist invariant gegen Anderungen der Amplitude 
erfafiter Signals. 



Die Anwendung der Erfindung fiir die Vorhersage oder Kontrolle 
eines Systems erfolgt derart, daii zunachst wie oben 
beschrieben aus der vergangenen Beobachtung und der Kenntnis 
■ der aktuellen Moden der IST-Zustand des Systems erfaiJt wird, 
der ggf . eine Mischung entsprechend dem Ergebnis der 
Driftsegmentierung darstellt. Der - IST-Zustand entspricht einem 
dynamischen System f. Die Vorhersage bedeutet, daii das System 
f auf den momentanen Zustand x angewendet wird und sich daraus 
die Vorhersage fur den unmittelbar folgenden Zustand y ergibt. 
Die Kontrolle bedeutet, daJi aus dem IST-Zustand die Abweichung 
von einem SOLL-Zustand ermittelt und aus der Abweichung eine 
geeignete Regelstrategie abgeleitet wird. 

Der Vorteil der Vorhersage und Kontrolle besteht darin, daB in 
komplexen Systemen (z. B. bei der Erfassung chemischer 
Reaktionen in einem Reaktor) , die ggf. die Messung lediglich 
weniger MeBgrofien erlauben, die fiir sich allein aufgrund von 
Mehrdeutigkeiten oder wegen systemimmanenten Verzogerungen 
keine direkten Riickschlusse auf den Systemzustand und ggf. 
bestehende gemischte Zustande erlauben, dennoch detaillierte 
Inf ormationen uber das System abgeleitet werden konnen. So 
, ^ _ ;ann bei dem Beispiel mit einer chemischen Reaktion aus der 
■^^c^^ierf indungsgemaBen Erfassung bespielsweise der makroskopischer 
thermodynamischer Zustandsgroflen eine optimale "Regelstrategie 
abgeleitet werden, die die Dosierung bestimmter 
Reaktionspartner umfaflt. 
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Patent ansprUche 

1 . Verf ahren zur Erf assung der Moden eines dynamischen 
Systems mit einer Vielzahl von Moden s^, die jeweils einen Satz 
a{t) charakteristischer Systemparameter besitzen, wobei eine 
Zeitreihe mindestens einer Systemvariablen x(t) einer 
Mcdellieriing unterzogen wird, die dazu eingerichtet ist, in 
jedem Zeitabschnitt einer vorbestimmten Mindestlange filr jede 
Systemvariable x(t) ein vorbestimmtes Vorhersagemodell fur 
2ine entsprechende Systemmode s^ zu erfassen, 

dadurch gekennzeichnet , dafi 

nach der Modellierung der Zeitreihe eine Drif tsegmentierung 
erfolgt, bei der in jedem Zeitabschnitt, in dem das System von 
einer ersten Systemmode Si zu einer zweiten Systemmode Sj 
ubergeht, eine Folge von gemischten Vorhersagemodellen g^ 
erfaiit wird, die durch eine lineare, paarweise Oberlagerung 
der Vorhersagemodelle f^ ^ der zwei Systemmoden s^^^ gegeben ist. 

2. Verfahren gemafi Anspruch 1, bei dem die Modellierung eine 
Schaltsegmentierung ist . 

Verfahren gemali Anspruch 2, bei dem die 
Schaltsegmentierung durch Simulierung einer 

Trainings-Zeitreihe des Systems oder der zu untersuchenden 
Zeitreihe mit mehreren, untereinander in Konkurrenz stehenden 
Vorhersagemodellen erfolgt , 

4. Verfahren gemaB Anspruch 3, bei dem die Vorhersagemodelle 
durch jieuronale Netzwerke oder andere Modelle zur Schatzung 
von Funktionen gebildet werden, die jeweils fiir eine Mode s 
charakteristisch sind und nach vorbestimmten Trainingsregeln 
um die Beschreibung der einzelnen Elemente der Zeitreihe 
konkurrieren . 



5 Verfahren gemaB einem der Anspriiche 1 bis 4, bei dem die 
Folge gemischter Systemmoden g^ aus den Vorhersagemodelle f^^^ 
und Interpolationsparametern a, b gemaiJ 
gi = a(s)fi(s)(x) + b (s) f-j(3) (x) ermittelt wird. 

6. Verfahren gemali Anspruch 5, bei dem die 

Interpolationsparameter gemaii 0 < a(s) < 1 und b(s) = 1 - a{s) 
gewahlt sind. 

7. Verfahren gemaB Anspruch 6, bei dem die Werte a(s) auf 
eine bestimmte Auf losungszahl R beschrankt und/oder 
aquidistant sind. 

8. Verfahren gemaB einem der " vorhergehenden Anspriiche, bei 
dem die Folge gemischter Vorhersagemodelle g^ erfaBt wird, 
indem fiir jeden Zeitschritt mit jedem der moglichen 
Vorhersagemodelle jeweils eine Vorhersage ermittelt wird, 
woraus sich eine zeitabhangige Vorhersagematrix ergibt, aus 
der ein mittlerer Vorhersagef ehler fur willkiirlich gewahlte 
Segmentierungen ableitbar ist, wobei die gesuchte Folge 
gemischter Vorhersagemodelle g^ die Segmentierung mit dem 
geringsten Vorhersagef ehler bzw. der maximalen 
Wahrscheinlichkeit ist . 

Verfahren gemaB Anspruch 8, bei dem die Suche nach der 
Segmentierung mit dem geringsten Vorhersagef ehler- mit einer 
dynamische Programmiertechnik erfolgt, die aquivalent zum 
Viterbi-Algorithmus fur Hidden-Markov-Modelle ist, wobei eine 
optimale Sequenz von Vorhersagemodellen unter Verwendung einer 
minimierten Cost-Funktion C* der Vorhersage ermittelt wird und 
die Segmentierung induktiv aus der Sequenz von 
Vorhersagemodelle abgeleitet wird. 

10. Verfahren gemaB einem der vorhergehenden Anspruche, bei 
dem nach der Drif tsegmentierung ein zusatzlicher Schritt der 



Reduzierung der Zahl der zur Modellierung verwendeten 
Vorhersagemodelle erfolgt, bei dem die Zahl der 
Vorhersagemodelle sequentiell jeweils verbunden mit einer 
Bestimmung des mittleren Vorhersagef ehlers so weit reduziert 
wird, dafl eine weitere Verringerung der Zahl der 
Vorhersagemodelle eine Erhohung des Vorhersagef ehlers 
bedeutet. 

11. Verfahren gemaB -einem der- -vorhergehenden Anspriiche, bei 
dem die Zeitreihe mindestens einer der Systemvariablen x(t) 
eine Zeitreihe von physiologischen Parametern umfaflt, die mit 
der Mackey-Glass-Verzogerungs-Dif f erentialgleichung 
^ = -O.lx(t) + ±111111^ 

beschrieben wird. 



12. Verfahren gemaB einem der Anspriiche 1 bis 11, bei dem die 
Zeitreihe mindestens einer der Systemvariablen x(t) eine 
Zeitreihe von physiologischen Parametern umfafit, die fur. den 
Ablauf von Schlaf- und Wachmoden charakteristisch sind. 



13. Verfahren gemafi Anspruch 12, bei dem die physiologischen 
Parameter EEG-Signale umfassen. 

14. Verfahren gemaB einem der Anspriiche 1 bis 10, bei dem die 
^.Zeitreihe mindestens einer der Systemvariablen x(t) eine 

'"^S:.-;^^2eitreihe von Sprachsignalen umfafit. 



ct3 




O 

o 



o 
o 

CsJ 
CN 



O 
O 
O 
CVJ 



o 
o 

00 



o 
o 

CO 



o 
o 



CO lO CO 



CO lO CO C\J T- 



o 
o 

CsJ 



bb 



